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' Abstract. We have assembled a sample of 37 radio-loud quasars that have been imaged with the Hubble Space Telescope in 

order to investigate their black hole masses, accretion rates, and the structure of their accretion disks. The black hole masses 
\l . were estimated from the luminosities of the host galaxies, and the accretion powers were extrapolated from the emission-line 

' luminosities. The majority of the quasars have masses in the range Mbh ~ 10* — 10^ A/q . Their accretion rates, M ~ 0.01 — 1 

times the Eddington rate, suggest that most of the objects possess standard optically thick, geometrically thin accretion disks, 
in some cases perhaps accompanied by an optically thin advection-dominated component. The coexistence of strong radio 
' emission and a standard disk conflicts with recent models for jet formation. We discuss modifications of the standard model 

, that can resolve this discrepancy. We find there is a strong correlation between the accretion rate and the extended radio 

QQ ■ luminosity. This lends support to the idea that the extended radio emission is somehow linked to the accretion disk. Lastly, we 

' combine the present sample of radio-loud quasars with the sample of BL Lac objects studied by Wang, Staubert, & Ho (2002) to 

CO . reevaluate the unification picture for radio-loud active galactic nuclei. Consistent with current ideas for the unification of radio- 

' loud sources, we find that flat-spectrum radio quasars and FR II radio galaxies indeed seem to belong to the same population, 

I as do BL Lac objects and FR I radio galaxies on the opposite end of the luminosity spectrum. However, some members of the 

low frequency-peaked BL Lac objects may be more closely associated with FR II rather than FR I radio galaxies. We describe 
I ' how the various subclasses of radio-loud sources can be viewed as a continuous sequence of varying accretion rate. 

o ■ 

5— ( ' 

■ Ksy words. BL Lacertae objects: general — galaxies: active — galaxies: jets — galaxies: nuclei — galaxies: Seyfert — 

^ ' quasars: general 
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1 . Introduction lobes. On smaller scales, Celotti & Fabian (1993) and Celotti, 
■ Padovani, & Ghisellini (1997) used radio data from very long- 
ed; According to the radio-loudness criterion proposed by basehne interferometry to estimate the kinetic luminosities in a 
Kellermann et al. (1989), quasars whose ratio of radio (5 GHz) j^j-gg sample of radio-loud objects. Both attempts reached the 
to optical (B band) flux exceeds 10 comprise a special pop- ^^me conclusion: the relativistic jet is linked with the accre- 
ulation of radio-loud objects (hereafter RLQs). As reviewed ^^^^ jisk. A more direct connection between the jet and the 
by Urry & Padovani (1995), the powerful radio emission is accretion disk has be found in the RLQ 3C 120 (Marscher et 
thought to originate from relativistic jets expeUed by the accre- 2002). However, how a relativistic jet actually forms from 
tion disk around a central supermassive black hole (BH). One j^e accretion disk remains an outstanding puzzle (Meier 2001 ; 
of the unsolved problems in astrophysics is the exact physical Blandford 2002) 
relation between the jet and the accretion disk. Two observa- 
tional approaches have been taken to investigate this relation The situation in Galactic BH candidates suggests that the 
in RLQs and active galactic nuclei (AGNs). On large (kpc to formation of the jet is somehow linked with flie accretion rate. 
Mpc) scales the relation between the jet radio/kinetic luminos- The radio emission from these objects becomes quenched when 
ity and the naiTow-line luminosity has been examined by many the system is in a soft (high accretion rate) state, as seen in 
aufliors (e.g., Baum & Heckman 1989a, 1989b; RawHngs & GX 339-4 (Fender et al. 1999), GRS 1915H-105 (Harmon et 
Saunders 1991). It has been found that the narrow-line lumi- al. 1997; Fender et al. 1999), GRO J16655-40 (Harmon et al. 
nosity is correlated with the kinetic luminosity from the radio 1995), Cynus X-3 (McCollough et al. 1999), and Cygnus X-1 
(Brocksopp et al. 1999). These objects may be instructive for 
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Table 1. The Sample of Radio-Loud Quasars 



No. 


Object 


Name 


Type 


z 


Mb (host) 


Mk(QSO) 


Reference 




Reference 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 


0133+207 


3C47 


S 


0.425 


-23.57 


-24.27 


2 


26.71 


6 


2 


0137+012 


PHL 1093 


S 


0.258 


-23.74 


-23.67 


2 


27.07 


7 


3 


0137+3309 


3C 48 


S 


0.367 


-25.03 


-25.92 


2 


27.52 


6 


4 


0202-765 


PKS 


S 


0.389 


-23.32 


-25.18 


2 






5 


0312-770 


PKS 


F 


0.223 


-23.98 


-24.60 


2 






6 


0340+048 


3C 93 


S 


0.357 


-23.93 


-23.71 


2 


25.56 


8 


7 


0405-123 


OF- 109 


F 


0.573 


-23.97 


-27.87 


1 


27.16 


6 


g 


0454-22 


PKS 


S 


0.533 


-21.39 


-26.79 


1 


27.57 


9 


9 


0736+01 


OI+061 


F 


0.191 


-23.78 


-24.23 


2 


24.32 


6 


10 


0837-120 


3C206 


S 


0.198 


-23.29 


-24.23 


2 


26.03 


6 


11 


0903+169 


3C215 


S 


0.412 


-23.28 


-24.28 


2 


26.51 


6 


12 


1004+130 


4C 13.41 


s 


0.240 


-24.26 


-25.82 


2 


26.15 


6 


13 


1020-103 


OL-133 


s 


0.197 


-23.36 


-23.47 


3 


24.83 


13 


14 


1138+0003 


LBQS 


7 


0.500 


-23.97 


-24.35 


4 






15 


1217+023 


ON+()29 


F 


0.239 


-23.71 


-24.40 


3 


25.45 


6 


16 


1218+1734 


LBQS 


S" 


0.445 


-23.46 


-23.74 


2 


25.93 


14 


17 


1222+1235 


LBQS 


s" 


0.412 


-24.02 


-24.24 


2 


25.96 


14 


18 


1226+023 


3C273 


F 


0.158 


-24.44 


-27.39 


2 


27.04 


6 


19 


1230-0015 


LBQS 


s 


0.470 


-23.14 


-24.64 


4 






20 


1250+568 


3C 277.1 


s 


0.321 


-22.96 


-23.29 


2 


26.73 


6 


21 


1302-102 


OP- 106 


s 


0.286 


-24.34 


-26.13 


2 


25.70 


6 


22 


1309+355 


Ton 1565 


F 


0.184 


-23.82 


-24.73 


2 


24.14 


6 


23 


1425+267 


Ton 202 


s 


0.366 


-23.65 


-25.69 


2 


25.78 


6 


24 


1512+37 


4C 37.43 


s 


0.371 


-23.86 


-25.58 


2 


26.28 


6 


25 


1545+210 


3C 323.1 


s 


0.266 


-23.53 


-25.14 


2 


26.42 


6 


26 


1548+ 114A 


4C+11.50 


F 


0.436 


-22.40 


-23.86 


2 


25.06 


10 


27 


1641+399 


3C 345 


F 


0.593 


-23.75 


-26.35 


1 


26.48 


11 


28 


1704+608 


3C 351 


s 


0.371 


-24.79 


-26.16 


2 


26.76 


6 


29 


2128-123 


OX- 148 


F 


0.501 


-22.82 


-27.11 


5 






30 


2135-147 


PHL 1657 


S 


0.200 


-23.43 


-24.16 


2 


26.34 


6 


31 


2141 + 175 


OX 169 


F 


0.213 


-23.38 


-24.79 


2 


26.33 


6 


32 


2201+315 


4C 31.63 


F 


0.295 


-24.70 


-25.98 


2 


26.54 


6 


33 


2247+140 


4C 14.82 


s 


0.237 


-23.52 


-24.10 


2 






34 


2348+0210 


LBQS 


? 


0.504 


-23.75 


-25.13 


4 






35 


2349 -014 


PB 5564 


s 


0.173 


-24.27 


-24.02 


2 






36 


2351-0036 


LBQS 


F 


0.460 


-23.02 


-24.18 


4 






37 


2355 -082 


PHL 6113 


s 


0.210 


-23.62 


-23.06 


3 


28.01 


12 



Notes — Col. (1) mimber; col. (2): object; col. (3): common name; col. (4): type (based on the spectral index between 1.4 and 5 
GHz; col. (5): redshift; col. (6) i?-band absolute magnitude of the host galaxy; col. (7): ij-band absolute magnitude of the QSO; 
col. (8): references for M/? ; col. (9): extended radio power (in units of W Hz~^) at 5 GHz; col. (10): references for i^r,f*Hz- 
REFERE.NCIiS— (1) Boyce, Disney, & Bleaken 1999; (2) Hamilton, Casertano, & Turnshek 2002; (3) Dunlop et al.''2()03; (4) 
Hooper, Impey, & Foltz 1997; (5) Boyce et al. 1998; (6) Xu, Livio, & Baum 1999; (7) Rector, Stocke, & EUingson 1995; (8) 
Bogers et al. 1994; (9) Aldcroft, Elvis, & Bechtold 1993; (10) Hutchings et al. 1996; (11) White 1992; (12) Lister, Gower, & 
Hutchings 1994; (13) Brotherton 1996. (14) Hooper, Impey, Foltz, & Hewett 1996 

based on the index between 0.365 and 4.85GHz (Gregory & Condon 1991 and Douglas et al. 1996) 

based on the index between 4.85 and 8.4 GHz (Gregory & Condon 1991 and Ref 14). 



lation between the "states" of the accretion disk and jet pro- 
duction in supermassive BH accretion systems is very poorly 
understood. 

There are two impediments in solving this problem. First, 
the mass of the central BH is highly uncertain. BH masses de- 
rived from fitting the "big blue bump" (Shields 1978; Malkan 
& Sargent 1982), for example, are coupled to the accretion rate, 
which is not known independently. Second, the observed con- 
tinuum of RLQs, even at optical and ultraviolet wavelengths, 
may be significantly contaminated by nonthermal emission 
from the relativistic jet. With the exception of objects with 
prominent big blue bumps, it is nontrivial to estimate the ther- 
mal component of the continuum. Thus, it is difficult to re- 
liably deduce the state of the accretion disk, which can be 
described conveniently by the dimensionless accretion rate 
m = M/MEdd, where M is the accretion rate, MEdd = 
2.2 X 10"** (t//0.1)"^ (Mbh/Mq) Mq yr^^ is the Eddington 
accretion rate defined by MEdd = -^Edd/'VC'^j -^Edd = l-26x 
10^^ (MBH/Af©) erg is the Eddington luminosity, and r\ is 
the radiative efficiency (Frank, King, & Raine 1992). 

The situation has improved dramatically with the availabil- 
ity of Hubble Space Telescope (HST) observations of quasar 



host galaxies. Relatively reliable luminosities for the host 
galaxies are now available. From studies of nearby inactive 
galaxies, it has been established that the mass of the BH scales 
roughly linearly with the luminosity (or mass) of the bulge 
component of the host galaxy (Kormendy & Richstone 1995; 
Magorrian et al. 1998; Ho 1999; Kormendy & Gebhardt 2001). 
Thus, we can use the luminosities of the quasar host galaxies, 
which are largely dominated by the bulge component, to esti- 
mate BH masses for quasars. 

Accretion disks have complicate structures for different m 
and the viscosity a (Chen et al. 1995). For a disk with a low 
accretion rate, three possible solutions have been suggested: 
optically thin advection-dominated accretion flows (ADAF; 
Narayan & Yi 1994), advection-dominated inflow-outflows 
(ADIOs; Blandford & Begelman 1999) and convection- 
dominated accretion flows (CDAF; Narayan, Igumenshchev & 
Abramowicz 2000; Quataert & Gruzinov 2000). Whether low 
m necessarily leads to low r], however, remains controversial. 
When < m < 1, the disk has a standard optically thick, 
geometrically thin structure (Shakura & Sunyaev 1973, here- 
after SS), which has been extensively studied. When m > 1, 
an opticaUy thick ADAF or super-Eddington accretion flow 
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Fig. 1. Distribution of redshift, Mi{(host), and Mii(QSO) for the present sample. 
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(SEA) with slim geometry is formed (Begelman & Meier 1982; 
Abramowicz et al. 1988; Chen & Taam 1993). Such a disk has 
two possible structures: an inhomogeneous structure (split into 
rings) due to the photon-bubble instability (Begelman 2002), 
or a homogeneous disk with photon trapping (Wang & Zhou 
1999; Ohsuga et al. 2002). Once the mass of the BH and the 
accretion luminosity are known, we can constrain the state of 
the accretion disk, which provides possible clues to jet forma- 
tion. 

Recently, a number of authors have suggested that flat- 
spectrum radio quasars, based on their very high kinetic lu- 
minosities, are powered by super-Eddington accretion flows 
(GhiselUni & Celotti 2001; Cavahere & D'Elia 2002; Maraschi 
& Tavecchio 2003). It would be highly desirable to test this and 
other hypotheses by obtaining a reliable estimate of the accre- 
tion rates of RLQs. This is the purpose of this paper. 

2. The Sample 

The most reliable measurement of quasar host galaxy luminosi- 
ties currently come from HST observations. Our analysis uses 
RLQs with available //5r imaging, based mainly on three sam- 
ples (see Table 1): (1) Hamilton et al. (2002), (2) McLure et 
al. (2001) and Dunlop et al. (2003), and (3) Boyce et al. (1998, 
1999). The B and V magnitudes are converted into the R band 
following Fukugita, Shimasaku, & Ichikawa (1995). We adopt 
the results from Hamilton et al. (2002) if the objects have been 
measured several times in order to decrease the systematic er- 
rors. Unfortunately, systematic errors in host galaxy measure- 
ments are difficult to quantify, but uncertainties on the order 
of 0.1-0.3 mag are not unexpected. Table 1 also lists the abso- 
lute magnitudes of the quasar nuclei, based on spatial decom- 
position of the HST images. Finally, we give published values 
of the extended radio luminosity. Following standard practice, 
the objects are divided into two subclasses according to tiieir 
radio spectral index n, defined by oc v^: flat-spectrum ra- 
dio quasars (FSRQs) are those with n > —0.5, whereas steep- 



spectrum radio quasars (SSRQs) have n < —0.5, where n is 
generally calculated between 1.4 and 5 GHz. There are two 
objects without spectral index, but they do not change our con- 
clusions (see Table 2 and 3). 

The redshift range of the sample is 0.16 < z < 0.57, and 
the host galaxy absolute magnitudes span —21.4 > Mb. > 
—25.0 (Fig. 1). For our present application, we use the empir- 
ical BH mass-bulge luminosity relation of McLure & Dunlop 
(2002): 



log (Mbh/Mq) = -0.5Mfl - 2.96. 



(1) 



For ease of comparison with literature results, distance- 
dependent quantities assume a Hubble constant of Hq = 50 km 
s~^ Mpc~^ and a deceleration parameter of qo = 0.5. 

3. States of the Accretion Disk 

The observed optical-UV continuum of RLQs may be con- 
taminated by nonthermal emission from the relativistic jet, 
rendering it difficult to estimate the thermal contribution to 
the continuum, and hence to estimate accretion rates. On the 
other hand, quasars, independent of their radio strength, usu- 
ally have strong, well-observed broad emission lines, which 
leads to the possibility of using the broad-line region (BLR) 
emission-line strength to indirectly estimate the ionizing lumi- 
nosity. Following Celotti et al. (1997), this is the approach we 
wiU take here. Celotti et al. (1997) showed fliat the BLR lumi- 
nosities in RLQs are roughly equal to the kinetic luminosities 
of their jets. We assume that the BLR clouds are most likely 
photoionized by the thermal emission from the accretion disk, 
rather than by nonthermal emission from the jet. This assump- 
tion is supported by that the beamed radiation from the jet has 
great difficulties to ionize the clouds distributed within a large 
solid angle to produce a reasonable range of ionization param- 
eter for emission lines. There is also ample evidence that BLR 
emission lines are photoionized by a thermal continuum (see, 
e.g., Netzer 1990). The underlying assumption here should be 
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Table 2. Line Luminosities 



No. 


Lines 


Reference 


log Flixx (erg s ^ cm 


log Luminosity (erg s ^ ) 


(1) 


(2) 


(3) 


(4) 


(5) 


1 


C IV, Lya 


15 


— 11.74 


45.22 


2 


Hq, H/3, H7 


1 


— 11.42 


45.09 


3 


H7, C IV, Mg 11, Lya 


16 


— 11.49 


45.34 


4 


H/3,H7, Mg II 


5 


-11.93 


44.95 


5 


C IV, Lya 


7,11 


— 11.41 


44.96 


6 


Ha 


6 


— 12.60 


44.21 


7 


Ha, H/3, C IV, Mg II, Lya 


7,9,10,11,12 


— 11.21 


46.03 


8 


Hp, C IV, Lya 


9,11 


— 11.64 


45.53 


9 


Ha, H/3, H7, C IV, Lya 


8,11 


-11.80 


44.43 


10 


Ha, H/3, H7, C IV, Lya 


8,10,11 


— 11.64 


44.63 


1 1 


H7, Mg II 


16 


— 11.95 


44.98 


12 


Ha, H/3, C IV, Mg 11, Lya 


8,10,11 


— 11.67 


44.77 


13 


Ha, H/3, H7 


13 


— 11.79 


44.47 


15 


Ha, H/3, H7, C IV, Lya 


8,11,7 


— 11.87 


44.57 


18 


H/3, H7, C IV, Mg II, Lya 


3,14 


— 10.28 


45.79 


20 


T T TT /O TT Tx r "X jff ^ Tx 

Ha, H/3, H7, C IV, Mg II 


8,15 


— 11.96 


44.75 


21 


H/3, C IV, Lya 


7,9,11 


— 11.51 


45.09 


22 


C IV, Lya 


11 


— 11.50 


44.70 


23 


Ha, H/3, C IV, Mg II, Lya 


10,11 


— 11.66 


45.16 


24 


Ha, H/3, C IV, Mg II, Lya 


7,9,10,11 


— 11.59 


45.25 


25 


H/3, C IV, Mg II, Lya 


7,9,15 


— 11.11 


45.42 


26 


H/3, H7 


2 


— 12.34 


44.64 


LI 


rip, ±17, IV, i^ya 


Q 1 1 
o,i I 






28 


Ha, H/3, C IV, Mg II, Lya 


7,9,10,11 


-11.63 


45.21 


29 


H/3, H7, C IV, Mg II, Lya 


2,7,11,12 


-12.04 


45.08 


30 


H/3, C IV, Lya 


7,9,11 


-11.40 


44.87 


31 


H/3, C IV, Lya 


7,9,11 


-11.69 


44.64 


32 


Ha, H/3, H7, C IV, Lya 


7,8,9,10,11 


-11.31 


45.31 


33 


H/3 


4 


-12.22 


44.21 


35 


H/3 


4 


-11.25 


44.89 



Notes — Col ( 1 ) number; col (2) emission lines; col (3) references; col (4) the derived total line fluxes; and col (5) the 

derived total line luminosity. 

REFERENCES— (1) Baldwin 1975; (2) Bergeron & Kunth 1984; (3) Corbin 1992; (4) Corbin 1997; (5) Danziger & 
Goss 1983; (6) Eracleous & Halpern 1994; (7) Gondhalekar 1990; (8) Jackson & Browne 1991; (9) Marziani et al. 
1996; (10) Neugebauer et al. 1979; (11) Osmer, Porter, & Green 1994; (12) Oke, Shields, & Korycanski 1984; (13) 
Scarpa & Falomo 1997; (14) Stockton & MacKenty 1987; (15) Wills et al. 1995; (16) Wilkes et al. 1999. 

quite safe in statistical test. 

Table 2 lists emission-line fluxes we were able to assem- 
ble for our sample. We choose prominent broad emission lines, 
such as Hq, H/3, H7, C iv, Mg 11 and Lya, that account for a 
large fraction of the total line luminosity. The overall similarity 
of the emission-line spectra of quasars suggests that the under- 
lying physical conditions of their BLR are quite similar (Netzer 
1990; Francis et al. 1991; Boroson & Green 1992; Zheng et 
al. 1997). To first order, therefore, the composite spectrum of 
quasars (Francis et al. 1991) may be used as a reasonable tem- 
plate to calculate the total BLR luminosity (Celotti et al. 1997). 
Then, from the relative contribution of individual lines to the 
total BLR luminosity, we obtain the total line luminosity I/unes 
(Table 2). 

The fraction ^ of the thermal disk emission reprocessed by 
the BLR clouds is roughly equal to the cloud covering factor. 
This assumption should be statistically valid for large samples 
since in steady state the energy absorbed by the BLR clouds is 
approximately that radiated as Une emission. Thus, the thermal 
emission from the accretion disk is 

-^disk = C lines- (2) 

According to Netzer (1990), ^ k, 0.1. 

We recognize that, for any individual object, variability 
and internal reddening may introduce uncertainties into iunes- 
Source-to-source variations in the covering factor or reprocess- 
ing efficiency may also add additional scatter. It is hoped, how- 
ever, that the overall statistics of a large sample wiU not be 
severely affected and that they are at least internally consistent. 



Following Wang, Staubert, & Ho (2002), we define the line 
accretion rate and its dimensionless form as 

A lines ines 
— 2~' ^ ~ "r • 

(r i-Edd 

The relation between A and the dimensionless accretion rate m, 
for an optically thin ADAF, is given by (Wang et al. 2002) 

m = 2.17 X lO-'ao.sCi'^'Ai^^ (4) 

This expression assumes that the total luminosity from the 
ADAF is Xdisk oc a-^MsHW^ (Mahadevan 1997), ao.3 = 
a/0.3, = ^/O.l, and A_4 = A/10~^. A necessary condi- 
tion for the presence of an optically thin ADAF is m < a^ 
(Narayan, Mahadevan, & Quataert 1998). Equation (4) can 
then be rewritten as 

Al = 1.72 X 10^3^_ia§.3. (5) 

Optically thin ADAFs require A < Ai. 

Accretion disks in the optically thick, geometrically thin SS 

regime obey 

m=^ = 10CiA, (6) 
and the condition 1 > m > a^ gives 

A2 = 9.0 X 10-3^_ia2 3, (7) 

above which a standard disk can exist. It is interesting to note 
that in the transition region between Ai and A2 the accretion 
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Table 3. Results 



Object 



10 




Log X 

Fig. 2. The distribution of (a) A and (h) A. The solid histograms plot 
the FSRQs, and the dotted-dashed histograms show the SSRQs. The 
dotted histograms show the total distribution for all RLQs (FSRQs 
and SSRQs). The distribution of A is divided into four regions, corre- 
sponding to different states of the accretion disks: (1) A < Ai (pure 
ADAF); (2) Ai < A < A2 (ADAF+SS); (3) A2 < A < A3 (SS); (4) 
A > A3 (SEA). 



flow may be in a hybrid state in which a standard disk coexists 
with an ADAF. The possibility of hybrid states in AGN accre- 
tion disks has been discussed by, among others, Quataert et al. 
(1999), Rozanska & Czemy (2000), Ho et al. (2000), and Ho 
(2002a). A transition from an SS disk to an ADAF is possible 
(Gu & Lu 2000), perhaps via evaporation (Liu et al. 1999). The 
transition radius depends on the accretion rate, the BH mass, 
and viscosity. The disk structure, however, is comphcated in 
such a regime, largely due to uncertainties in the viscosity. 
When the accretion rate reaches to > 1, we have 



A3 = O.lC-i, 

in which case we use 



Er = 



(8) 



(9) 



log Mbh 

(Me) 



log A 



log X Disk 





8.82 


24.27 


-1.71 


SS 


2 


8.91 


24.14 


-1.92 


SS 


3 


9.56 


24.39 


-2.32 


SS-fADAF 


4 


8.70 


24.00 


-1.85 


SS 


5 


9.03 


24.01 


-2.17 


SS-^ADAF 


6 


9.01 


23.26 


-2.90 


ADAF 


7 


9.02 


25.08 


-1.10 


SS 


8 


7.73 


24.58 


-0.31 


SEA 


9 


8.93 


23.48 


-2.60 


SS+ADAF 


10 


8.69 


23.68 


-2.16 


SS+ADAF 


11 


8.68 


24.03 


-1.80 


SS 


12 


9.17 


23.82 


-2.50 


SS-^ADAF 


13 


8.72 


23.52 


-2.35 


SS-^ADAF 


15 


8.89 


23.62 


-2.43 


SS-^ADAF 


18 


9.26 


24.84 


-1.57 


ss 


20 


8.52 


23.80 


-1.87 


SS 


21 


9.21 


24.14 


-2.22 


SS+ADAF 


22 


8.95 


23.75 


-2.35 


SS+ADAF 


23 


8.86 


24.21 


-1.81 


ss 


24 


8.97 


24.30 


-1.82 


ss 


25 


8.81 


24.47 


-1.49 


ss 


26 


8.24 


23.69 


-1.70 


ss 


27 


8.91 


24.49 


-1.58 


ss 


28 


9.44 


24.26 


-2.33 


SS-^ADAF 


29 


8.45 


24.13 


-1.47 


ss 


30 


8.76 


23.92 


-1.99 


ss 


31 


8.73 


23.69 


-2.19 


SS+ADAF 


32 


9.39 


24.36 


-2.18 


SS+ADAF 


33 


8.80 


23.26 


-2.69 


SS+ADAF 


35 


9.18 


23.94 


-2.39 


SS+ADAF 



to describe Ihe accretion state, since rh > in the super- 
Eddington regime in which photon trapping efficiently lowers 
the radiated luminosity (Wang & Zhou 1999; Ohsuga et al. 
2002). 

Our estimation of the parameter A are subject to at least 
three sources of uncertainty. First, the total fine luminosity is 
estimated from only a limited number of emission lines. Since 
the relative strengths of the main strong lines are relatively well 
determined, and quasars generally have fairly uniform spectra 
(e.g., Francis et al. 1991; Zheng et al. 1997), this should not 
introduce a large source of uncertainty. More worrisome is our 
assumption of a constant cloud covering factor (10%), which 
in reality is Ukely to have an intrinsic dispersion. And third, 
even if we assume that the host galaxy magnitudes are per- 
fectly measured, the Mji — Mbh relation still has a scatter of 
~ 0.4 dex (McLure & Dunlop 2002). For concreteness, we as- 
sume that any individual value of A may have an uncertainty 
of ~ 0.5 dex. Although this is not negUgible, fortunately it is 
not sufficiently large to obscure gross statistical trends in our 
sample. However, we caution against a literal interpretation of 
the individual values of A and A listed in Table 3. 

The three critical values of A define four regimes in accre- 
tion states, as shown in Figure 2b. The average values of A 
for FSRQs and SSRQs are, respectively, (log A>fsrq = 24.0 
and (log A)ssRQ = 24.1. The corresponding values for A are 
(logA)FSRQ = -2.20 and (logA)ssRQ = -2.05. Given the 
uncertainties just discussed, the two subsamples appear to be 
statistically indistinguishable. 

Wang et al. (2002) present the distributions of A and A for 
BL Lac objects and conclude that all BL Lacs have ADAFs. 
Comparing those objects with the RLQs presented here, we 
find that the distributions of A overlap significantly, increasing 
smoothly from BL Lacs to RLQs. This was noted by Scarpa 
& Falomo (1997), who used only Mg 11 as the luminosity in- 
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dicator. Our result considerably strengthens theirs. A similar 
behavior is seen in the distribution of A. As we discuss later, 
these trends suggest that an evolutionary sequence, character- 
ized by different accretion states, may connect BL Lac objects 
to RLQs. 

The normalized mass accretion rates, m = (,~^X « lOA, 
span over a relatively large range, from 10"-^ to lO^ *^. Some 
of the FSRQs in our sample fall in the range of objects stud- 
ied by Cavaliere & D'Elia (2002), but we suggest that some of 
FSRQs, like SSRQs, may have hybrid ADAFh-SS disks. With 
the similarity in m between FSRQs and SSRQs, the only dif- 
ference between FSRQs and SSRQs is geometrical orientation, 
which has been suggested (Urry & Padovani 1995). 

Cavaliere & D'Elia (2002) suggest that an evolutionary se- 
quence may be defined for blazars, depending on their dimen- 
sionless accretion rates. FSRQs have m w 10^^, whereas BL 
Lac objects have m w 10~^, the two populations separated by 
a significant gap. By contrast, our analysis suggests that BL 
Lacs and RLQs define a continuous sequence in accretion rate. 
Although our conclusion is based on two heterogeneous sam- 
ples, we note that it is in general agreement with the evolution- 
ary scenario proposed by Bottcher & Dermer (2002). 

The states of the accretion disk listed in Table 3 can be 
tested by Fe Ka line observations. Differences in the X-ray 
properties between radio-loud and radio-quite quasars have 
been discussed by Sambruna, Eracleous, & Mushotzky (1999) 
and Reeves & Turner (2000). The most important findings in 
these studies is that the Fe Ka line tends to have complex pro- 
files. Only a few RLQs have broad Fe Ka profiles, most be- 
ing undetected (Sambruna et al. 1999; Reeves & Turner 2000). 
There is also a clear trend that more radio-loud objects tend to 
have lower Fe Ka equivalent widths. By analogy with the in- 
terpretation for the relativistically broadened Fe lines in Seyfert 
nuclei (e.g., Nandra et al. 1997), RLQs with strong broad Fe 
Ka lines may be in the SS state, while those lacking such a 
feature may be in a pure ADAF or in an ADAFh-SS state, ac- 
cording to the classification suggested in this study. This in- 
terpretation, although appealing, may not be unique. Recall 
that among Seyfert galaxies and radio-quiet quasars there ex- 
ists an inverse correlation between Fe line strength and X-ray 
luminosity (Iwasawa & Taniguchi 1993), an effect that can 
be attributed to luminosity-dependent changes in the ioniza- 
tion state of the surface of the disk (Reeves & Turner 2000). 
Alternatively, perhaps strong relativistic jets physically modify 
the structure of the innermost part of the disk, thereby altering 
or suppressing the iron-line emitting region (see § 4.1.3). 

We note that values of Ai and A2 are sensitive to a. The 
most likely value of a is 0.3 (Narayan, Mahadevan & Quataert 
1998). However if a = 0.1, it will follow that Ai = 1.91 x 
10-''^_iag.i and A2 = 1.0 x IQ-^^^iag.i. All the disks Hsted 
as hybrid models of SS and ADAF will be SS disk in Table 3. 
In any case, this does not affect the conclusion that most of the 
objects have the standard accretion disks. 

There is only one possible super-Eddington accretor (PKS 
0454—22) in the present sample. This nicely agrees with the 
result from an independent test by Wang (2003), in which 
super-Eddington accretors have been searched in the present 
sample and no super-Eddington accretor has been found. He 



used the theoretical spectra emergent from slim disks (Wang 
et al. 1999) and the empirical reverberation relation to find a 
limit relation bewteen the BH mass and H/3 width as Mbh = 
(2.9- 12.6) xIO^Mq (i;FWHM/103kms-i)® '^^inAGNswith 
super-Eddington rates. With the help of equation (1), the BH 
masses can be obtained, he found all the objects in the present 
sample are located below this limit. 

It is worth noting that, with the exception of 3C 273, our 
sample lacks very powerful FSRQs. This is an observational 
selection effect. It is difficult to measure the host galaxies of 
the most powerful quasars because of their very bright nuclei. 
We will later return to the issue of whether there exists super- 
Eddington RLQs that appear as very powerful FSRQs. 

4. Discussion 

4. 1. Jet Formation from the Disk 

Two processes have been widely advocated for jet forma- 
tion in AGNs: (1) energy extraction from the spin of the BH 
(Blandford & Znajek 1977, hereafter BZ) and (2) energy ex- 
traction from a disk wind (Blandford & Payne 1982, hereafter 
BP). Several variants of the BZ and BP models have been pro- 
posed. The results of the present study allow us to set some 
constraints on the theoretical models. In their investigation of 
the BZ process in the regime of the SS disk, Livio, Ogilvie, & 
Pringle (1999) conclude that both the BZ and BP power, even 
for Kerr BHs, are negligible compared with the radiative output 
from the accretion disk itself. The reason is that geometrically 
thin disks cannot produce a strong enough poloidal magnetic 
field for the energy extraction to be effective. 

4.1 .1 . Model of Ghosh & Abramowicz 

Ghosh & Abramowicz (1997) assume that the poloidal mag- 
netic field can be parameterized as Bp oc aPmax, where Pmax 
is the maximum pressure in the inner region if the accre- 
tion rate rh > 10~^(a_2rn9)~^/^, with a_2 = a/10~^ and 
mg = Mbh/IO^Mq. In such a region, the radiation pressure 
dominates over the gas pressure and is independent of the ac- 
cretion rate. Then, the BZ power is given by 

Lbz = 2 X lO^'^mg/ crgs-\ (10) 

where j — JeH/CGMgjjC^^) is the specific angular momen- 
tum. We can see that the BZ power is independent of the ac- 
cretion rate, but very sensitive to the spin. The extended ra- 
dio power is expected to be proportional to Lbz (Meier 1999). 
For our sample. Figure 3 shows that the extended radio pow- 
ers span ~ 5 orders of magnitude, whereas the BH masses 
are confined to within roughly a factor of 10. This means that 
the BHs should have spins in the range 10"^ < j < 1-0. 
Such a large range of values for the spin, however, is not ex- 
pected in RLQs (Wilson & Colbert 1995). The model of Ghosh 
& Abramowicz (1997), therefore, seems to be disfavored by 
the present data. Additionally, the results of Lacy et al. (2001) 
from FIRST Bright Quasar Survey show that the radio emis- 
sion depends on the accretion rate (oc m^ °), which does not 
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support Equation (10). As pointed out by Livio et al. (1999), 
Equation (10) overestimates the BZ power because the poloidal 
magnetic field should be Bp w {H/R)B^ rather than simply 

4.1 .2. Model of Meier 

Meier (2001) considers energy extraction from the spin of 
the BH and from the disk itself, for both an SS disk and 
an ADAF. His model uses the middle-region solution for the 
SS disk and a self-similar solution for the ADAF; it consid- 
ers both Schwarzschild and Kerr BHs. In this model the jet 
power for a Schwarzschild BH is only Ljet = 10^^ '' erg s~^ 
if Mbh = IO^Mq and to = 0.1, whereas Ljet = lO^^ '^ erg 
s"'^ for a Kerr BH with the same mass and accretion rate. The 
present data clearly do not support this model because many 
RLQs have SS disks that support much higher luminosities. 
When the accretion disk is in an ADAF state, for example, an 
accretion rate of m = 10~^-^ onto a Schwarzschild BH ra- 
diates ~ 10"'^^ erg s~^, which may generate a radio power 
of 10^^ '^e_2W Hz^^ at 1 GHz in the extended lobes, where 
e_2 = e/10~^ is the fraction of the total jet power radiated 
in the radio band (Meier 1999). This can account for the radio 
emission from BL Lac objects shown in Figure 3; it also im- 
plies that the BHs of BL Lac objects may not necessarily have 
large spins. The jet power from an ADAF surrounding a Kerr 
BH is approximately lO^'''^ erg s"'^, corresponding to a radio 
power of 10^^-^e_2W Hz~^ at 1 GHz. This can account for 
RLQs in the ADAFh-SS regime. But how do we explain those 
RLQs that have much higher radio powers and that contain pure 
SS disks? 



4.1 .3. An Improved Model 

From the above discussion, it is evident that the power output of 
the jet, whether generated through the BZ or BP processes, de- 
pends both on the model adopted for the accretion disk and the 
parameterization used for the poloidal magnetic field. Below 
we derive an improved formula for the jet power. 

Following Meier (2001), the magnetohydrodynamic 
(MHD) power of the jet is given by 



1 



i'jet — ^7:^B„Rr,fl, 



2 n4ri2 



32c 



(11) 



where Bp is the poloidal magnetic field and Oq is the Keplerian 
angular velocity at the radius -Ro> where the jet develops. Both 
the BZ and the BP process need a strong poloidal magnetic 
field Bp, which is produced by the toroidal magnetic field B^. 
Livio et al. (1999) showed that 



Bp^l- 



Ba 



(12) 



where H is the height of the disk at radius R and the index n ~ 
1. Since H/R < 1 for an SS disk, it is generally believed that 
Bp <C Bfj,, and hence that the extracted energy must be small 
(Livio et al. 1999; Meier 1999, 2001). However we should note 
that H/R « 0.01 is only valid in the middle and outer regions 



of the disk; in the inner region, H/R « 1. Meier (2001) argues 
that this region is thermally unstable. He only uses the middle- 
region solution in his calculations, and consequently obtains a 
reduced power from the BZ process. 

The thermal instability arises from insufficient cooUng. It 
is well known that such an instability in the innermost regions 
of the disk can be suppressed by advection, which acts as an 
efficient cooling process (Abramowicz et al. 1988; Chen & 
Taam 1993). Advection mainly converts the heat of the accret- 
ing gas into entropy, rather than radiation. Consequently, the 
sound speed (Cg) and the scale height H are increased. It would 
be of interest to include the role of advection cooling in the 
calculation of the BZ power. 

We follow Chen's (1995) treatment of advection cooUng in 
determining the structure of the disk. The disk is assumed to 
undergo Keplerian rotation at an angular velocity Qk, and we 
assume conservation of mass, conservation of angular momen- 
tum, and vertical equilibrium. Advection is included as local 
cooling, and we neglect the corona. The energy balance be- 
tween the heat production (Qvis) and coohng can be written 
as 



^vis — Qrad "I" ^adv) 



(13) 



where Qrad and Qadv represent radiative and advection cool- 
ing, respectively. The viscous heating rate per unit area is 
Qvis = 3Mf2^/5(/47r, where the correction terms / ~ 1 and 
5 w 1 are due to the inner boundary condition and the pseudo- 
Newtonian potential (Chen 1995). Advection cooling can be 
expressed as Qadv = Mc'^C/2ttR-^, ( = 5.5 being the ad- 
vection factor. With the diffusion approximation and the as- 
sumption of radiation pressure dominance, we have Qrad = 
4acT^/3Kes5^ = ^c^kCs/kss, where S is the surface density 
of the disk, a = 7.56 x 10~^^ is the radiation density constant 
and Kcs = 0.34 is the electron-scattering opacity. With the help 
of the equations for mass and angular momentum conservation 
and vertical equihbrium. Equation (13) reduces to 



2cf]K 3M 

Cs - -;— = 0. 

Kes 47r 



(14) 



The first term represents advection cooling, which is only im- 
portant at small radii or for high accretion rates, and the second 
represents radiative coohng. Neglecting the first term yields the 
standard solution for an SS disk, whereas neglecting the sec- 
ond gives the solution for an optically thick ADAF (or SEA). 
Advection coohng helps to suppress thermal instabilities that 
develop in the irmer regions of the disk. 

The solutions of the disk give the sound speed 



as = 0.09 mljr°-^ jyTTsSrhV" 



1 



where Us = Cg/c and to_i = rh/O.l, the pressure 
P = 1.76 X lO'^ a'^rh-iTOg ^a^V"^ dyne cm"^, 



and the vertical half-thickness 
H 
'R 



(15) 



(16) 



(17) 
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Fig. 3. The distribution of extended radio power at 5 GHz versus the (a) A, (b) A, and (c) Mbh- The filled circles show the FSRQs, and SSRQs 

are denoted as open stars. The open and semi-filled circles represent high frequency-peaked and low frequency-peaked BL Lacs, respectively; 
the radio data for these objects are taken from Urry et al. (2000), and values for A, A, and Mbh are taken from Wang et al. (2002). The diagonal 
line in panel (c) is the empirical division between FR I and PR II radio galaxies in the diagram of extended radio power at 1.4 GHz versus 
optical R-betad absolute magnitude, as presented by Ledlow & Owen (1996); see text for details. 

As shown by Xu et al. (1999), the mass outflow rate from 

the jet is given by 

/ TT\ 2<5-3 

Mjet « (^-^j M, (21) 

where 6 is a constant from 1.7-3.4. For Kerr BHs, the char- 
acteristic radius tq = 1.5 leads to « 0.65r^"-^, and 
hence H/R 0.65 in the jet-forming region. We then have 
A-fjct ~ (0.19 - 0.84)M: the mass outflow rate of the jet is a 
large fraction of the accretion rate. If the jet-formation region is 
an ADAF, the outflow rate will be even higher than this because 
H/R « 1. We note that such a large outflow rate may have 
implications for observations of the broad Fe Ka line in radio- 
loud AGNs. The region of the disk that emits the red wing of 
the broad line roughly coincides with the region where the rel- 
ativistic jet is formed. If strong winds accompany jets, they can 
substantially alter the inner structure of the disk. This may be 
at least partly responsible for the observed diversity of Fe-line 
strengths and the complexity of their profiles in RLQs (Reeves 
& Turner 2000). 

We note that Equation (21) is not self-consistently obtained 
from the disk model. Such a high outflow rate may take place 
in ADIOS models. Here we stress that the Fe Ka line in RLQs 
may be strongly affected by the outflow. The case of 3C 120 
adds credibility to the picture that jet formation influences the 
innermost region of disk: "dips in X-ray emission are fol- 
lowed by ejections of bright superluminal knots in the radio 
jet" (Marscher et al. 2002). The detectabiUty of the Fe Ka line 
in radio-loud objects may also be influenced by contamination 
of the X-ray continuum by non-thermal (beamed) emission. 

4.2. Empirical Evidence for a Disk- Jet Connection 

The correlation between the extended radio and line emis- 
sion has been explored by a number of authors (e.g., Baum & 



where r = R/Rg and Rg = GMbu/(^. With flie help of 
B'^/Stt f» aP and Equation (12), we have 

Bp = 2.02 X loWj^rrig^^'^al/'^r-'^ Gauss. (18) 

Following Meier (2001), we take the characteristic ra- 
dius -Ro = 7Rg and Rq = l.bRg for a Schwarzschild 
and Kerr BH, respectively, and the angular velocity of the 
metric in the Boyer-Lindquist frame ui = —g4,t/(j<b<b ~ 
0.3j (GMbh/c^) ^- We finally get the jet power from the ac- 
cretion system. For a Schwarzschild BH, 

Ljet = 8.60 X 10^^ mgrh-i erg s"\ (19) 

and for a Kerr BH, 

Ljet = 1.87 X lO''^ TO9TO_i(l + l.lj + 0.3j2) crgs-^ (20) 

Considering that the bolometric luminosity from accretion is 
1.38 X lO^^mgrh-i erg s^^, the ratio of the jet power to ac- 
cretion power for a Schwarzschild BH is roughly 4 x 10~^, 
whereas for a Kerr BH it is ~1.0. 

The output powers given by Equations (19) and (20) are 
different from those found by Ghosh & Abramowicz (1997) 
and Meier (2001). First, the energy output based on Equations 
(19) and (20) depends on the accretion rate, whereas the model 
of Ghosh & Abramowicz (1997) does not. Second, the output 
is much higher than that given by Meier (2001) for the standard 
accretion disk model. 

The MHD power for a Schwarzschild BH is due to the BP 
process. For a Kerr BH, the BP process also depends strongly 
on the spin because the metric itself drags the MHD fluid, 
which then contributes to the rotation of the magnetic field; 
thus. Equation (20) covers both the BP and BZ processes. The 
above levels of MHD power, which is extracted from the BH 
and the disk in the standard SS regime, are sufficient to account 
for the highest radio luminosities seen in RLQs. 
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Heckman 1989a, 1989b; Rawlings & Saunders 1991; Rawlings 
1992; Falcke, Malkan, & Biermann 1995; Xu et al. 1999; Cao 
& Jiang 1999; Willott et al. 1999). The present sample allows 
us not only to further test this correlation, but presents us with 
an opportunity to explain it in the context of our improved 
model. Figures 3a and 3b show that the extended radio emis- 
sion is evidently quite closely linked with A and A. The linear 
regression for the RLQs alone is 

log Lg^Hz = (28.05 ± 0.78) + (0.95 ± 0.38) log A, (22) 

with a Pearson's correlation coefficient r = 0.56, significant 
at the level of p = 3.8 x 10~^ according to Student's t test. 
Similarly, 

log^sGHz = (-8.65 ± 11.03) + (1.45 ±0.46) log A, (23) 

with r = 0.71 and p ^ 6.2 x lO'^. 

For comparison, we include the BL Lac objects studied 
by Wang et al. (2002) that have extended radio luminosities 
given by Urry et al. (2000); there are 12 high frequency-peaked 
(HBL) and 15 low frequency-peaked (LBL) objects. The BL 
Lac objects smoothly extend the sequence defined by the RLQs 
toward lower values of Xsghz, A, and A. The correlations are 
now enhanced, becoming 

log^IcHz = (27.57 ± 0.29) + (0.68 ± 0.10) log A, (24) 

with r = 0.65 and p = 2.6 x 10"^°, and 

logilcHz = (4.18 ±2.70) + (0.92 ±0.11) log A, (25) 

with r = 0.76 and p = 2.1 x 10~^^. In the above correlation 
analysis of L^^^ versus A, we excluded the upper limits for 
the BL Lac objects. A survival analysis including these upper 
limits shows log LI^Hz = (27.60±0.30) + (0.64±0.11) logA 
with Spearman's p = 0.667. These correlations are striking in 
that they provide fairly direct evidence that the extended jet 
emission is somehow closely coupled to the central engine via 
the accretion disk. The most direct interpretation is that the ex- 
tended radio emission is fed by relativistic jets that are tightly 
linked to the accretion disk (Rees 1984). 

Following Meier (2001), we assume that 1% of the MHD 
energy will be radiated as radio emission with a spectral index 
of —1. Then, for a Kerr BH and using the relation A = rj^M, 
L^^tjj^ = 3.74 X 1026e_2m9m_i W Hz"!, or 

L^gjj^ = 3.74 X lO^^e-amgA W Hz'^ 

= 2.69 X lO^e-aA WHz'^ (26) 

These expressions are in rough agreement with the observed 
correlations. 

4.3. A Power Limit for Radio-Loud Quasars? 

In our sample, only PKS 0454—22 is a super-Eddington ac- 
cretor {En = 10° ^^). The most powerful RLQs, which have 
powers in excess of 10*^ erg s~^, potentially may be super- 
Eddington accretors (Ghisellini & Celotti 2001). In this section 
we consider the upper luminosity limit of RLQs as a threshold 
set by super-Eddington accretion. 



9 

For a BH accreting well above the Eddington rate, much of 
the dissipated energy will be advected into BH because most 
of the photons are trapped in the flow by Thomson scattering 
(Wang & Zhou 1999; Ohsuga et al. 2002). Photon trapping ef- 
ficiently controls the radiated luminosity; the radiated luminos- 
ity from the disk depends very weakly on the accretion rate, 
^disk oc log M, tending to be saturated. For a self-similar disk 
(Wang & Zhou 1999), Equation (15) reduces to 



as ~ 0.65r-°•^ (27) 
and the MHD power is given by 

Ljet = 8.05 X 10*^ mimg erg s"^ (28) 
for a Schwarzschild BH and 

Ljet = 3.76 X 10^^ rhimg erg s"^ (29) 



for a Kerr BH, where nii = m/lQ. The power is still propor- 
tional to the accretion rate, although the radiation from disk 
itself is saturated. In the super-Eddington regime, the differ- 
ence in radio power between a Kerr and Schwarzschild BH is 
smaller. 

From the total rotational energy due to the spin, -Egpin ~ 
1.6 X lO^^mgj^ erg, it follows that the spin-down timescale 
due to the BZ process is 

Tspin-down = 1.9 X 10^ iTi^^ yr. (30) 

This is much shorter than a typical merger timescale of 10^ 
yr (Wilson & Colbert 1995), if the spin angular momentum 
comes from galaxy collisions. If the powerful jet is fed by 
a Schwarzschild BH via super-Eddington accretion, the e- 
folding timescale of the BH growth will be 

Tgj^^ = Mbh/M = 4.0 X lOVf Vr- (31) 

The lifetime of the powerful jet from an SEA is rather short un- 
less there is continuous super-Eddington fueling. Thus, super- 
Eddington RLQs ought to be very rare. 

4.4. AGN Unification 

A useful diagnostic for radio galaxies, as first shown by Owen 
& Ledlow (1994; see also Ledlow & Owen 1996 and Owen, 
Ledlow, & Keel 1996), is the plot of extended radio luminos- 
ity versus optical absolute magnitude. A diagonal line well 
separates the two radio populations: at a given optical lumi- 
nosity (which in the case of radio galaxies is dominated by 
the host galaxy), FR lis have much higher radio powers than 
FR Is (Fanaroff & Riley 1974). The physical interpretation of 
the Owen-Ledlow line has been much debated. Bicknell (1995) 
and Gopal- Krishna & Witta (2001) suggest that it may arise 
naturally from the interaction of the jet with the interstellar or 
intracluster medium. Or perhaps the two classes reflect intrin- 
sic differences in their central engines, such as BH spin (Baum, 
Zirbel, & O'Dea 1995; Meier 1999) or state of the accretion 
disk. Ghisellini & Celotti (2001) recently suggested that the di- 
viding fine corresponds to a transition in the accretion mode, 
from an SS disk (FR II) to an opticaUy thin ADAF (FR 1). 
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Urry et al. (2000) also examined BL Lac objects and RLQs 
in the Mr — i^Sio diagram, but the data for their RLQs are 
mostly not from observations. Figure 3c gives a new rep- 
resentation of the Owen-Ledlow diagram for RLQs and BL 
Lacs, where we have substituted the optical luminosity on the 
abscissa with BH mass A/bh- We have plotted our sample of 
RLQs along with the sample of BL Lac objects studied by Urry 
et al. (2000). The diagonal line is the Owen & Ledlow (1994) 
line for radio galaxies, adapted from Ledlow & Owen (1996) 
after converting to our distance scale, the radio luminosities 
from 1.4 GHz to 5 GHz assuming a spectral index of —1, and 
Mfj to Mbh using Equation (1). 

Several features are noteworthy in this diagram. First, the 
Owen-Ledlow Une that divides FR I and FR n radio galax- 
ies seems to cleanly bisect the population of beamed radio- 
loud objects. Nearly all HBLs fall below the line, and all but 
two FSRQ lies above the Une. Second, the population of LBLs 
straddle the divide, occupying the region between HBLs and 
FSRQs. Third, as with the FSRQs, essentially all SSRQs hne 
above the line. And fourth, at least some BL Lac objects, 
mostly of the LBL variety, appear to be more closely associ- 
ated with FR 11 radio galaxies; this is in agreement with the 
suggestion of KoUgaard et al. (1996) that some radio- selected 
BL Lac objects may be beamed FR II sources. 

The availability of accretion rates enables us to consider 
radio-loud AGNs in a more unified, physical context. Wang 
et al. (2002) show that most BL Lac objects have accretion 
rates below the critical value to w that demarcates the 
transition from a standard SS disk to an optically thin ADAF. 
Moreover, they find that HBLs have characteristically lower ac- 
cretion rates than LBLs. This study demonstrates that nearly all 
RLQs have to > a^. It is important to note that, at least for 
the sample considered here, most of the sources occupy a rela- 
tively narrow range of BH masses, ~ 10^-^ Mq with a spread 
of ~0.5 dex. The various subpopulations of radio-loud AGNs 
can therefore be cast into a physical, perhaps evolutionary, se- 
quence described by a single variable — accretion rate. For the 
beamed population, to decreases along the sequence FSRQ 
LBL HBL. When viewed at large angles to our line of sight, 
the sequence becomes FR II (FR 11 or FR 1) FR 1. Since 
the structure of the accretion flow is governed by to, the varia- 
tion in TO along the sequence directly translates into a variation 
in the dominant mode of accretion. The most powerful RLQs 
and FR II radio galaxies may contain super-Eddington disks or 
optically thick ADAFs. The majority of RLQs and FR lis have 
SS disks. The low-power HBLs and some LBLs have pure op- 
tically thin ADAFs. Objects that occupy the fuzzy boundary 
near the Owen-Ledlow line, which comprise some LBLs and a 
minority of FR lis, may have a hydrid structure consisting of 
an SS disk plus an optically thin ADAF. 

The basic elements of the above evolutionary scenario, 
based on complementary but different lines of evidence, have 
already been advocated by a number of authors (e.g., Ghiselhni 
et al. 1998; GhiseUini & Celotti 2001, 2002; GhiselUni, Celotti, 
Costamante 2002; Bottcher & Dermer 2002; CavaUere & 
D'Elia 2002). Bottcher & Dermer (2002) and Cavaliere & 
D'Elia (2002) argue that RLQs evolve into BL Lacs rather than 
BL Lacs into RLQs (Bottcher & Dermer 2002). If the evolu- 



tionary direction is from RLQs to BL Lacs, the BH masses 
of RLQS should be systematically smaller than those of BL 
Lacs. However, the significant overlap in the BH mass dis- 
tributions of the two groups (Fig. 3c; see also Falomo 2003) 
does not seem to support this single evolutionary direction. 
Given the current uncertainties in BH mass determinations for 
AGNs, it is probably premature to draw strong conclusions. 
Here, we merely wish to note that both evolutionary paths can, 
and probably do, occur. The detection of quasar remnants in 
nearby galaxies in the form of inactive massive BHs conclu- 
sively demonstrates that the cosmic evolution of AGN activity 
must be episodic (Richstone et al. 1998). The duty cycle for ac- 
cretion is short, such that any individual massive BH is likely to 
have been activated (and deactivated) many times since it was 
formed. 

5. Conclusions 

We attempt to constrain the accretion rates and the mode of ac- 
cretion in RLQs using BH masses estimated from HST imag- 
ing of their host galaxies and accretion luminosities derived 
from their emission-line spectra. The quasars have large BH 
masses, A'/bh ~ 10® — lO'^ Mq, and substantial accretion 
rates, Af « 0.01 — 1 times the Eddington rate. Their accretion 
rates suggest that most of the objects possess standard optically 
thick, geometrically thin accretion disks, perhaps in combina- 
tion with an optically thin ADAF in some cases. Our sample 
has a deficit of objects with super-Eddington disks, although 
this is likely due to a selection effect in the sample. 

Our data can also be used to test current ideas on jet for- 
mation. We present an improved model for a standard optically 
thick, geometrically thin disk that incorporates the effects of 
advection cooling. Our model can explain the level of radio 
emission observed in RLQs. We discuss some consequences of 
our model for the interpretation of Fe Ka lines observed in the 
X-rays. 

Finally, we combine the present sample of RLQs with a 
recently analyzed sample of BL Lac objects to reevaluate the 
unification picture for radio-loud AGNs. The strong correlation 
found between the luminosity of the extended radio emission 
and the accretion rate supports the notion that radio jets are di- 
rectly coupled to the accretion disk. In agreement with prevail- 
ing ideas, our results support the proposition that FSRQs and 
BL Lacs are the beamed counterparts of FR 11 and FR 1 radio 
galaxies, respectively. However, we argue that some fraction of 
the low frequency-peaked BL Lac sources in fact may be more 
closely associated with FR II radio galaxies. The various sub- 
classes of radio-loud AGNs — ranging from high-luminosity 
sources (FSRQs and FR Us) on the one end to low-luminosity 
sources (BL Lacs and FR Is) on the other — can be viewed as 
a continuous sequence of varying accretion rate. 
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